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We present an analysis of exclusive charmless semileptonic B-meson decays based on 83 million 
BB pairs recorded with the BABAR detector at the T(45') resonance. Using isospin symmetry, we 
measure branching fractions B{B° ■k-(.+v) = (1.38 ± 0.10 ± 0.16 ± 0.08) x 10"" and B{B° 
p~ t^v) = (2.14±0.21±0.48±0.28) x 10""*, where the errors are statistical, experimental systematic, 
and due to form-factor shape uncertainties. We compare the measured distribution in , the 
momentum-transfer squared, with theoretical predictions for the form factors from lattice QCD 
and light-cone sum rules, and extract the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\Vub\ = (3.82 ± 0.14 ± 0.22 ± Q.lltoH) x 10"^ from B n£u, where the fourth error reflects the 
uncertainty of the form-factor normalization. 

PACS numbers: 13.20.He, 12.15.Hh, 12.38.Qk, 14.40.Nd 



The parameter \Vub\ is one of the smallest and least 
known elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) quark-mixing matrix A precise determina- 
tion of \Vub\ would significantly improve the constraints 
on the Unitarity Triangle and provide a stringent test 
of the Standard Model mechanism for CP violation. In 
this paper, we present a determination of \Vub\ from 
charmless semileptonic decays of B mesons with exclu- 
sively reconstructed final states, B huiv, where the 
hadronic state /i„ represents a tt^, tt^ , p^, or and £ 
represents e or p. Exclusive decays allow for kinematic 
constraints and more efficient background suppression 
compared to inclusive decays, but must rely on theoreti- 
cal form-factor predictions. Using isospin symmetry, we 
measure the branching fractions B{B'^ Tr~£'^v) |2]| 
and B{B'^ p~£^v) as a function of = (pg+p,^)'^, 
the momentum-transfer squared, and extract \Vub\ us- 
ing recent form-factor calculations based on light-cone 
sum r ulcs (LCSR) H 13 and unqucnchcd lattice QCD 
(LQCD) ^1^]. 

This measurement is based on a sample of 83 million 
BB pairs recorded with the BABAR detector 0] at the 
PEP-II asymmetric-energy e"'"e~ storage rings. The data 
correspond to an integrated luminosity of 75.6 fb^^ col- 
lected at the T{4S) resonance and 8.9 fb~^ recorded 40 
MeV below it. Simulated BB events are used to esti- 
mate signal efficiencies and shapes of signal and back- 
ground distributions. Charmless semileptonic decays 
are simulated as a mixture of three-body decays B — > 
Xuir^ {Xu = TT^rj,!]' , p,uj) based on the ISGW II quark 
model Decays to non- resonant hadronic states Xu 
with masses mx > SttItt are simulated following a pre- 
scription of Ref . 13 . 

We identify charmless semileptonic decays by a 
charged lepton with momentum \p^\ > 1.3 GeV |22|, a 
TT or p meson, and missing momentum |pmiss| > 0.7 GeV 
in the event. We identify p mesons via the decays 
p^ tt^tt" and p^ tt+tt^ with mass 0.65 < rriTrTT < 



0.85 GeV, rejecting candidates in which a charged track 
is identified as a kaon; both tt^ and p candidates are re- 
jected if a charged track is identified as a lepton. The 
charged lepton is combined with a 7r°, or tt^, p^ of 
opposite charge to form a "y" candidate; Y candidates 
are rejected if the lepton and an oppositely-charged track 
from the signal hadron are consistent with a. J/^ 
decay. 

The neutrino four- momentum, pi, = (S'missiPmiss), 
is inferred from the difference between the net four- 
momentum of the colliding-beam particles, Pbeams = 
(-E-boamsiPboams), and the sum of the four-momenta of 
all detected particles in the event. To reduce the ef- 
fect of losses due to the detector acceptance, we re- 
quire a total charge of the event of IQtotI < 1 and 
a polar angle of the missing momentum in the range 
0.6 < 6',niss < 2.9 rad. In addition, the missing mass mea- 
sured from the whole event should be compatible with 
zero. Because the missing-mass resolution varies linearly 
with the missing energy, we require \m'^^^^^/2E^iss\ < 
0.4 GeV. We compute the angle between the Y can- 
didate and the B meson, assuming zero missing mass, 
as cos6'by = (2^;* - M| - M^) / (2|p* |). Here 
MB,MYTEg,EY,PBiPY refer to the masses, energies, 
and momenta of the B and Y . Signal candidates are 
required to satisfy | cos^byI < 1.1, allowing for detector 
resolution and photon radiation. 

We restrict the momenta of leptons and hadrons in 
Y candidates to enhance the signal over backgrounds. 
For B ttIv, we require \p( \ + \Ph \ > 2.6 GeV; for 
B plv, 1.5|p;| -I- |p;*J > 4.2 GeV and \p^\ > 1.8 GeV. 
These criteria keep 99.8% (75%) of true B tt{p)Iv de- 
cays and reduce the B Xctv background by about 
10% (80%) after all other selection criteria. To suppress 
backgrounds from e+e" <M{<1 = u,d,s,c) and QED 
processes, we require at least five charged tracks in each 
event or, to increase the efficiency for i?+ 7r°£+zy, 
four tracks and at least two photons. We also require 
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L2 = X^i IpTIcos^^* < 1.5 GeV. Here the sum is over all 
tracks in the event excluding the Y candidate, and p* 
and 6* refer to the momenta and the angles measured 
with respect to the thrust axis of the Y . This require- 
ment removes over 95% of qq and 80% of i? — *■ Xciv 
background and retains about 50% of the signal in all 
modes. 

We discriminate against the remaining background us- 
ing the variables Ai? = {pb ' J'bcams — 

s/2)/Vi and 

TOES = ^/{s/2+PB ■ Pbeams)^/-E^oams ~ Pb ' whcrC ^ is 

the mass of the T{AS). Only candidates with |AiJ| < 
0.9 GcV and mEs > 5.095 GeV are retained. The to- 
tal signal selection efficiencies for the sum of electrons 
and muons are 3.5% and 2.4% for tt~£~^i/ and 7r°^+i^, 
0.53% and 1.1% for p-i+v and p°£+iy. [H We use a 
low-background sample of D*~t^v decays with 

D° i<r+7r~ or — *■ ii'+7r~7r° to compare the efficien- 
cies of each selection cut in data and simulation and find 
differences typically of a few percent. 

To extract the signal yields, we perform a binned ex- 
tended maximum-likelihood fit to the A£^ vs. toes dis- 
tributions of the four signal modes simultaneously. The 
fit takes into account statistical ffuctuations of both data 
and Monte-Carlo samples. We fit the relative propor- 
tions of the simulated signal and background samples to 
the data distributions in 5GeV^ or lOGeV^ intervals of 
q^. To improve the q^ resolution, we adjust \py\ so that 
AE = 0. The resulting resolution is small compared to 
the chosen intervals of q^ and can be described by the sum 
of two Gaussian functions of widths cri ~ 0.2 GeV^ (con- 
taining about 75% of signal events) and a2 — 0.5 GeV^. 

We use the isospin relations T{B'^ 7r~£+j/) = 
2T{B+ ^ TT^l+v) and T{B^ p-l+v) ^ 2T{B+ 
p^t^v) to reduce the number of fit parameters to nine: 
five for the signal yields in the five q^ intervals for 
B -kIv decays, three for the signal yields in the three 
intervals for B —f piv decays, plus one scale parameter, 
shared among all q^ intervals and signal modes, to fit the 
overall normalization of the B Xctv background. We 
classify signal candidates as "combinatoric signal" if the 
reconstructed lepton comes from the isospin-conjugate 
decay or the hadron is incorrectly selected. The fit uses 
common parameters for combinatoric signal and signal. 
The normalization of the simulated non-BB background 
is scaled separately for events with and p^ to match 
the off-resonance data. We smooth the distributions for 
this low-statistics background to reduce single-bin statis- 
tical fluctuations. 

Figures m and 121 show projections of the fitted A_B vs. 
Toes distributions for each q^ interval for B —> iriv and 
B piv, respectively. Integrated over the whole 
range, we observe 396 TT~i^v, 137 7r°£+j/, 95 p~i~^i^, 
and 98 p^i'^v decays. The resulting partial and total 
branching fractions are given in Tabled The fitted nor- 
malization of the B Xciv background is consistent 
with the measured total branching fraction The 



TABLE L Partial and total branching fractions B{B° — » 
Tv-e+u) (AB^) and B{B" p-£+u) (ABp) obtained from 
the simultaneous fit of the four signal modes, and signal ef- 
ficiencies, e,r and tp, averaged over charged and neutral B 
decays. The errors are statistical. 



Range 






Range 


ABp 




(GeV^) 


(10-*) 


(%) 


(GeV^) 


{W~^) 


(%) 


0-5 


0.30 ±0.05 


2.1 










0-10 


0.73 ±0.17 


0.70 


5-10 


0.32 ±0.05 


2.9 








10-15 


0.23 ±0.05 


3.8 


10-15 


0.82 ±0.10 


0.97 


15-20 


0.27 ±0.05 


3.5 








15-25 


0.59 ±0.07 


0.44 


20-25 


0.26 ±0.03 


3.3 








0-25 


1.38 ±0.10 


3.1 


0-25 


2.14 ±0.21 


0.72 



goodness-of-fit is evaluated using a x^-based comparison 
of the fitted AE vs. toes distributions and data, yielding 
X^/dof = 1.27. As a check, we have performed the fit for 
and p^ separately and obtain consistent results. 
The fit also allows us to study the q^ dependence of the 
form factors. In decays to pseudoscalar mesons there is 
only one form factor, (for low-mass leptons), and we 
can extract the shape of f+{q^) directly from the mea- 
sured q^ spectrum. We perform a fit with a function 
proposed by Becirevic and Kaidalov (BK) [lll |. 

/■ /„2^ CB{l-a) 

■'^^^ ' {l^q^/ml,){l-aq^/ml.) ' ^ ' 

where raB" = 5.32 GeV is the mass of the B* resonance, 
cb is a normalization factor, and a is a shape param- 
eter. Since we cannot measure the normalization, only 
a is meaningful. Leaving both cb and a free, we fit 
a = 0.61 ± 0.09, in agreement with LQCD results dill. 
For decays to vector mesons, there are three form fac- 
tors. The experimental uncertainties for B pt.u are 
still too large to measure these. Thus we have to rely on 
theoretical predictions. 

Figure O compares the distributions for B niv 
and B piv with the various form- factor calculations, 
which we implement by reweighting simulated signal 
events ^3 ■ We use probabilities to quantify the agree- 
ment: for B Triiy we obtain good agreement with the 
BK fit to the data, P{x^) = 35%; and the predictions 
of LCSRl 01, 38%; LQCDl 0, 14%; and LQCD2 H], 
35%; but only marginal agreement with the prediction 
of ISGW II 0, P(x^) < 1%. For B piu all calcula- 
tions 0, Q are compatible with the data within the large 
experimental uncertainties. 

The systematic errors in the extraction of the branch- 
ing fractions are listed in Table ^ The contributions 
from each q^ interval are conservatively treated as fully 
correlated and added linearly to obtain the uncertainty 
of the total branching fractions. Part of the q^ variation 
of the stated errors may be due to statistical variations 



6 




5.1 5.2 5.1 5.2 5.1 5.2 5.1 5.2 5.1 5.2 

m,3(GeV) 




-0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 -0.5 0.5 



A E (GeV) 

FIG. 1: (color online) Projected tties (a-e) and AE (f-j) distributions in five intervals of for the combined B — > Trliy modes. 
The projections are shown for signal bands —0.15 < AE < 0.25 GeV and rriES > 5.255 GeV, respectively. The error bars on the 
data points represent the statistical uncertainties. The histograms show simulated distributions for signal (white), combinatoric 
signal (white, dotted), cross feed from other B Xulv decays (hatched), B Xclv decays (light shaded/yellow) and non-BB 
background (dark shaded/blue). The normalizations of the signal and B — > X^iv background simulations have been scaled to 
the results of the maximum-likelihood fit. 



TABLE II: Relative systematic uncertainties of the partial and total branching fractions S(i?" n f.'^v) {AB-^) and B{B'^ 
t^u) {ABp) in the various bins. The total uncertainty in each column is the sum in quadrature of the listed contributions. 



Range (GeV^) 


0-5 


SABn/ABn (%) 
5-10 10-15 15-20 20-25 


0-25 


5ABp/ABp (%) 
0-10 10-15 15-25 0-25 


Track and Photon Reconstruction 


7.4 


5.7 


9.2 


3.5 


8.7 


6.8 


16.7 


10.8 


15.9 


14.2 


K% Production and Interactions 


8.8 


5.0 


7.3 


2.5 


4.8 


5.7 


12.7 


4.7 


10.9 


9.1 


Lepton Identification 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


2.3 


B — » Xc^v Background 


5.0 


3.9 


3.9 


4.3 


3.6 


4.2 


7.2 


1.8 


3.8 


4.2 


B Xuf-i' Background 


0.5 


1.5 


0.6 


2.2 


5.7 


2.1 


10.9 


9.1 


19.2 


12.5 


Non-B_B Background 


13.5 


2.4 


1.0 


2.2 


7.8 


5.6 


11.2 


0.9 


1.6 


4.6 


Nbb 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


B Lifetimes 


1.1 


1.4 


0.2 


0.5 


2.1 


1.1 


0.7 


0.6 


0.9 


0.7 


/+-//00 


0.7 


0.8 


0.3 


0.4 


1.0 


0.7 


1.6 


0.6 


0.4 


0.9 


Isospin Breaking 


0.1 


1.1 


1.8 


0.1 


0.1 


0.6 


6.4 


4.5 


0.9 


4.2 


Radiative Corrections 


0.8 


0.3 


0.2 


0.3 


1.5 


0.6 


0.5 


0.1 


0.9 


0.5 


Total Error 


18.7 


9.5 


12.8 


7.3 


14.8 


11.8 


28.1 


15.9 


27.7 


22.5 



in simulated samples. 

Uncertainties in the simulation of the reconstruction 
of charged particles and photons are evaluated by vary- 
ing the reconstruction efficiencies and the photon-energy 
resolution and are added in quadrature. In addition, 
most A'2 escape detection. The impact of A'° interac- 



tions in the calorimeter is estimated by varying in sim- 
ulation their detection efficiency and energy deposition. 
To assess the uncertainty of the K'^ production rate, we 
vary the inclusive branching fractions of K'^X, 
D° K^X, and D+ K°X within their published er- 
rors ^3 ■ these constitute the total uncertainty of the 
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-0.5 0.5 

A E (GeV) 



0.5 0.5 



FIG. 2: (color online) Projected m^s (a-c) and AE (d-f) 
distributions in three intervals of for the combined B — > plu 
modes. See the caption for Figure Q for details. 



0.3 



-a 



0.2 



0.1 



(a) 



ISGW II 
- LCSR 1 ■ 

- LQCD 1 
LQCD 2 

- BK Fit to Data 
_i_ 



10 15 
(GeV') 



0.6 



0.4 



0.2 







(b) 



-ISGW II 
- LCSR 2 
- Fit to Data 



20 25 



10 15 

(GeV^) 



20 25 



FIG. 3: (color online) Comparison of the differential decay 
rates as functions of for B — > irlu (a) and B — > plu (b) with 
various form-factor predictions. The data are background 
subtracted and corrected for efficiency and radiative effects. 
The error bars are statistical (inner) and statistical plus sys- 
tematic (outer). 



neutrino reconstruction, which is dominant. For lepton 
identification we use relative uncertainties of ±2% and 
±4% for electrons and niuons, respectively. 

The uncertainty of the B Xctv background is eval- 
uated by varying the B D/D*/D**£i/ branching frac- 
tions [131 and the B D*iv form factors [l3|- For 
the B — > Xu(.i^ background, we independently vary the 
branching fractions of B^ — > ujI'^v and B^ — *■ rjt^v 
within their published errors 0, 0| . We assume equal 
branching fractions for rjiv and rj' and use a relative 
uncertainty of 100% for the latter. We also vary the non- 
resonant contribution within the range allowed by the un- 
certainty of the total B X^tv branching fraction [lij . 



The impact of quark-hadron duality violation or weak 
annihilation effects have not been considered. We es- 
timate the uncertainty of the small remaining non-BB 
background by comparing simulation with off-resonance 
data and extract a normalization error of for elec- 

trons and ±25% for muons. 

The overall uncertainty of the number of produced B 
mesons is 1.1%. We take into account the uncertainties of 
the ratio of S lifetimes, tb± /tbo = 1.081±0.015 [l3|, the 

charged-to-neutral B production ratio /-I //oo = 1.044± 

0.050 01, and the potential effect of isospin breaking due 
to p° — Lj mixing [13 ■ We assign an uncertainty of 20% 
to the radiative corrections based on PHOTOS [l^ . 

The impact of the uncertainties of the B nlv 
form-factor shape on the measured branching fractions 
is negligible, whereas for the different B — > piv form- 
factor calculations we see variations of up to ±6% in 
B{B" n-e+iy) and ±13% in B{B° p-(+v). We take 
the full spread between calculations as the uncertainty of 
the dependence of the form factors. 



TABLE III: \Vub\ derived for B -> ntv and B -> plu sig- 
nal for various q^ regions and form-factor (FF) calculations: 
LCSRl 0, LQCDl II, LQCD2 g, LCSR2 0, ISGW II 0. 
For the cross feed from the other mode, we have used the BK 
fit to data for -kIv and LCSR2 for plv. Quoted errors are sta- 
tistical, experimental systematic, uncertainties of form-factor 
shape and form-factor normalization (no form-factor nor- 
malization uncertainties are available for plv). 







Range 


AC 


IKbl 






(GeV^) 


(ps-i) 


(10-^) 




TT FF 












LCSRl 





-15 


5.1±1.3 


3.27 ±0.16 ±0.19 ±0.10 


f 0.53 

-0.36 


LQCDl 


15 


-25 


1.5±0.4 


4.92 ±0.25 ±0.29 ±0.15 


fO.76 
-0.52 


LQCD2 


15 


-25 


2.0±0.5 


4.16 ±0.22 ±0.24 ±0.12 


4-0.72 
-0.47 


LCSRl 





-25 


7.7±2.3 


3.40 ± 0.13 ± 0.20 ± o.io1E;:^2 


LQCDl 





-25 


5.7±1.7 


4.00 ±0.14 ±0.23 ±0.12 


fO.78 
-0.49 


LQCD2 





-25 


6.1±2.1 


3.82 ± 0.14 ± 0.22 ± 0.1llo:g 


p FF 












LCSR2 





-15 


12.7 


2.82 ±0.18 ±0.30 ±0.18 




ISGW II 





-25 


14.2 


2.91 ±0.12 ±0.33 ±0.19 




LCSR2 





-25 


17.2 


2.85 ±0.14 ±0.32 ±0.19 





We extract \Vub\ (see Table from the partial 
branching fractions AS using the relation \ V,ih\ = 
^A6/(tboAC) , where tbo = (1.536 ± 0.014) ps [ij is 
the i?° lifetime and A^ denotes the predicted form-factor 
normalization in each interval. For q^ < 15 GeV^ we 
derive \Vub\ using LCSR calculations; for > 15 GeV^ 
we use unquenched LQCD. To extract \Vub\ from this 
measurement over the whole range, we extrapolate 
the LQCD results to low using the fits of the BK 
parametrization in Ref. [J, I3 and the LCSR results to 
high q^ using the parametrization given in Ref. 0. We 
adopt the uncertainties of the form-factor normalization 
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estimated in Refs. H S S H . 

In conclusion, we have measured the exchisivc branch- 
ing fractions B{B° ti"(+v) and B{B^ p-(+v) as 
a function of q^, and have extracted \Vub\ using recent 
form-factor calculations. We measure the total branch- 
ing fractions, 

B{B^ ^TT~tv) = (1.38 ± 0.10 ± 0.16 ± 0.08) X 10~^ 
B{B° ^ p-l+v) = (2.14 ± 0.21 ± 0.48 ± 0.28) X 10"'', 

where the errors are statistical (data and simulation), 
experimental systematic, and uncertainties of the form- 
factor shapes. As a consistency check, we have also 
measured the branching fractions for the charged and 
neutral niv samples separately, B{B^ ir^i'^v) = 
(1.41 ± 0.17 ± 0.17 ± 0.08) X 10-'', B{B+ Tr°e+v) ^ 
(0.70 ± 0.10 ± 0.08 ± 0.04) x lO"''. The ratio T{B° 
Tr-e+iy)/T{B+ Tr"e+v) = 2.21 ±0.41 is consistent with 
the assumed isospin relation within the quoted statistical 
uncertainty. 

For B tt£i^, the distribution agrees well with cal- 
culations based on LCSR and unquenched LQCD ^ 
0, but the data disfavor ISGW II 0. Instead of aver- 
aging results based on different calculations, we choose 
the measured form-factor shape and normalization of 
LQCD2 and quote 

iKbl = (3.82 ± 0.14 ± 0.22 ±0.1lt^:g) x 10"^ 

where the additional fourth error reflects the uncertainty 
of the form-factor normalization. The results arc con- 
sistent with previous measurements 0, ll^l, but have 
higher statistical accuracy, are less dependent on theo- 
retical form-factor predictions, and benefit from recent 
advances in theoretical calculations H> Hi - 
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